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Vacuum polarization in external fields is treated by way of calculating — exactly and then 
perturbatively — the phase of the quantum scattering matrix in the Shale-Stinespring approach 
to field theory. The link between the Shale-Stinespring method and the Epstein-Glaser renor- 
malization procedure is highlighted. 



1. Introduction 

Renormalization theory is bound to suffer a reappraisal in the light of the reconstruc- 
tion of Zimmermann's forest formula in Hopf-algebraic terms, together with the interpre- 
tation of the dimensional regularization method given by Connes and Kreimer [1] , and the 
birth of quantum field theory on noncommutative spaces [2, 3], together with the evidence 
that Yang-Mills theories on noncommutative manifolds are ultraviolet divergent — see [4] 
and the other references in that paper. 

Now, part of the advantages and new popularity (see, for instance, [5-8]) of the 
Epstein-Glaser renormalization method [ 9,10] stems from the fact that it is locally de- 
fined, and so in principle applies to models on nonflat manifolds. There is however some 
contention on whether, as claimed by some practitioners [11], the Epstein-Glaser method 
is a fundamental one. 

In cases like these, it sometimes helps to look at a simpler problem, for which an 
absolutely reliable method is known, to see how the marketed procedures fare in its respect. 

The chosen problem is that of vacuum polarization in external fields and the chosen 
reliable method is the Shale-Stinespring approach to linear quantum field theories. This 
is an entirely rigorous algebraic method; in [3] its prima facie applicability to "imple- 
mentable" theories on noncommutative as well as commutative spaces was brought to the 
fore. 

After reviewing the Shale-Stinespring theorem in Section 2, in the body of this article 
we show, by a refinement of its technique, that the phase of the scattering matrix is well 
defined and finite for implementable linear theories. We give explicit formulae for the 
phase. 

We then exploit in QED a perturbative version of this approach, which leads in a direct 
way to the formulas reached by Scharf [11] in his account of vacuum polarization in QED 
by the Epstein-Glaser method. Some simplification of his calculations results from using 
gauge-invariant variables. The contention by Scharf and followers that the Epstein-Glaser 
renormalization procedure is a fundamental one is vindicated to some extent. 

2. A reminder on Shale-Stinespring theory 

Let H = H + © 7i~ =: P+H © P-7i be a Hilbert space, graded by the projections P± 
on the positive and negative spectral subspaces for a free Dirac operator. Operators on TC 
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are presented in block form: 
We have in mind particularly the classical (or "first-quantized" ) scattering matrix 



S 



S ++ S+- 
S-+ S- 



i S++ \ „ / S+- 

PS '■= Oeven = I n c , qs := Oodd ~" 



Introduce a nomenclature for its even and odd parts: 

'++ \ a -« 

o S— ) ' qs — " oaQ " v S-+ o 

Similarly, 



s*_ si J' Ps, = l o si J' ?s, = lsi. 

Unitarity of S 1 gives the identities 

S , ++<S++ + S+-S+_ = S++S++ + S^_ + S-+ = P+, 
S—S j __ + S- + S f _ + = S t __S— + S j + _S+- = P_, 
S ++ Sl + + S+Sl_ = Sl + S+- + = 0, 

ssl_ + S-+sl + = sl_S- + + sl_s ++ = 0. 

It is clear that p^ 1 exists if and only if S++ and S are invertible as operators on 7i + 

and on 7i~ , respectively; this is the generic case, that always holds when S is close to the 
identity, and will be the only one considered in the sequel. We then define the skewadjoint 
operators 



Ts:=qsp ~ s ={s. + s-l "o "J' Ts:=Ts '-\-szls. + 

Consider the Fock space constructed on H with the new scalar product 

(v I <P) ■= (V+ I <P+) + (<P- \v-), (1) 

where rj± := P±rj. Let {4>k\ and {V>fc} denote arbitrary orthonormal bases for H + and TL~ , 
respectively; we shall abbreviate bk '■= b(4>k), dk '■= d(ifjk) in the notation of the "particle" 
and "antiparticle" annihilation operators, and similarly for the creation operators b\ : d\. 
For any operator A on Ti we have the quantum (or "second-quantized") counterpart, acting 
on Fock space: 

dA(A) := b j A ++ b + tfA + _d^ + dA_+b + :dA—d j : . 
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Here, for instance, :dA d^: denotes — ^ ■ k d\.{ipj \ A ipk)dj, the double colon meaning, 

as usual, a normally ordered product, and A^ d^ is ^Zj k b\{4>k I A + _ipj)dj; the other 

cases should be clear. 

This rule corresponds to the infinitesimal spin representation [12]; it is independent of 
the drafted ortho normal bases and makes sense only when A^ , A^ are Hilbert-Schmidt. 

The rule is mainly applied to selfadjoint operators, and yields (at least formally) 
selfadjoint operators in turn. For instance, in QED the free Dirac equation is written as 

d . n , _ d , 
where the Dirac matrices, say in the chiral representation, are given by 

To this classical selfadjoint operator D Q corresponds the quantum free Hamiltonian 

H := dA(D ) = b ] D b + :dD d ] : , 
which is a positive operator. The classical interaction Hamiltonian is of the form 

H{t) =e{A°(t)-a-A(t)), 

where e denotes the electromagnetic coupling constant and (A , A) =: A the electromag- 
netic vector potential, a real c-number function. Note the covariant form H(t) = e^p A, 
with A := / y fJ, A fJi . In the interaction picture one considers 

V(t) := e iDot H(t)e- iDot . (2) 

The quantum interaction Hamiltonian is then 

V(t) :=dA(V(t)) = tfV++(t)b + bW + _(t)d^ + dV_+(t)b + :dV—(t)d^: . (3a) 

This can be rewritten in terms of the formal fermion field as 

V(t) = J d 3 x-^(x)4(x)^(x): 1 (36) 

with x = (t, x) and the bar meaning the Dirac adjoint. For that, just write the fermion 
field in the form 

= Mk{x) + d\%l) k (x)), 
k 

where <j>k(x) = e~ lDot (f)k(x), and similarly for the ip : s. 
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For the quantum scattering matrix S we need instead the global spin representa- 
tion [12,13], which we call A. It is given by 



S := e ie A(S) = e ie |(0 in | out )| :expdA(J): = (0 in | out ) :expdA(J): . 



Here 0i n denotes the incoming vacuum, ou t := S 0i n , and 



( 



s+_s. 

szis_ + i-s: 



I := 



-l 



Again, this makes sense if and only if S-\ , S ^ are Hilbert-Schmidt, and then we say that 

S is implementable. The absolute value of the vacuum persistence amplitude (0i n | out ) is 
given by 



For our present purposes, this is the content of the Shale-Stinespring theorem [12, 14]. 
The phase 9 is in principle undetermined and conventionally taken equal to zero; this is 
all that is needed to compute transition probabilities. 

The phase of the vacuum persistence amplitude does matter physically, however: the 
current density is modified with respect to the free field situation by the vacuum polariza- 
tion effect (that bears on the radiative correction to the photon propagator in the nonlinear 
theory), and the interacting current density is found by functional derivation of S with 
respect to the gauge potential, in which the phase intervenes [15]. 

The question is then to find an appropriate and computable definition for the phase 
of the quantum scattering matrix. 

3. Computing the phase in the Shale— Stinespring framework 

The difficulty comes from the fact that the global spin representation is projective. 
Let U(s,t) be the classical unitary propagator in the interaction representation, which 
interpolates between the identity and S. Then U(s, t) solves the equation 



|(0 in | out )| = detail - T|) = dee/ 2 (S__Sl_) = det 1 / 2 (S ++ S{ + ) 
= det 1/2 (l - S + _Sl_) = det 1/2 (l - S_ + SL + ). 




(4) 



an explicit form being given by the Dyson expansion 




from which follow the propagator properties: 



U(t,t) = l, 



U(t,s)U(s,r) = U(t,r). 



4 



Here S = U(oo, — oo). 

Now, for unitary implementable operators Ui, U 2 , in general 

A(tfi)A(*7 2 ) = c(U u U 2 )A(U 1 U 2 ), 

where the cocycle c (with the vanishing phase convention) is given [12, 13] by 

c(Ui, U 2 ) = exp(zargdet 1/2 (l - T^fj/J) = exp^argdet 1 / 2 ^ 1 ^^^ 1 ))- (5) 

This is not a trivial cocycle because the determinants of the pu operators are not individu- 
ally defined in general. Assume that the interpolating family U(s, t) is implementable and 
(strongly) differentiable with respect to its parameters — this will happen if the external 
field is sufficiently well behaved. Then 

A(U(s, t))A(U(t, r)) = c(s, t, r)A(U(s, r)), (6) 

with an obvious notation. Still, c(t,t,r) = c(s,t,t) = 1. On the other hand, it can be 
shown [12] that 

A(U(s,t)) = V(t). 

98 s=t 

We now seek to redefine A(U(s,t)) by multiplying a phase factor e ie ^ s,t " > so that the 
new quantum family U(s,t) := e i9 ^ s ^ A(U(s, t)) fulfils 

U(t,t) = l; V(s,t)V(t,r) = V(s,r), (7) 

just like the classical propagator. If we manage that, then 6>(+oo, —00) will have every right 
to be called the phase of the quantum scattering operator S. Let c(s, t, r) =: exp(i£(s, t, r)). 
Differentiating equation (6) , one gets: 



V(t)A(U(t,r)) = i^- s 



c(s,t,r)A(U(t, r)) + i%-A(U(t,r)). 

s=t 01 



Then, we redefine 



U( S ,£):=exp^ A ^(\,r,t) dr^A(U(s,t)), 



which clearly satisfies 

U(s,t) = l-i I V(u)V(u,s)du. 



This equation is the quantized version of (4) and sports the same kind of iteration solution: 

U(s,t) = l + J2(-*) n V(ti) / V(t 2 )..-/ V(t n )dt n ...dt 2 dt 1 . (8) 
1 Jt Jt Jt 



n=l 



In the present case, the quantum Dyson expansion is rigorous: although V is an unbounded 
operator, it is a pretty tame one. Let E m denote the projector on states containing at 
most m particles. Then, in view of (3), Vi? m is a bounded operator from E m into E m+ 2, 
and the norms of 



U(s, t) :n E r , 



V(ti) / V(t 2 ) 



V(t n ) dt n ... dt 2 dt x E r , 



can be easily estimated. To see that, one introduces the norm 

III^WIII := ||Kven(t)|| + ||Kdd(*)||2, 

where the latter is the Hilbert-Schmidt norm. By continuity and uniform boundedness, 
|||^(^)||| < o(s, r) holds for some finite function a(s, r), when r < t < s, and is not difficult 
to check that there are constants C m such that 

II V(t)E m \\ <C m a{s,r). 

Consult [16, 17] for precise analyses of these bounds (encompassing also the boson case). 
On the other hand, from the integral equation (4), 

|||t7(s,r)||| < l + {s-r)a(s,r). 

Putting both inequalities together, one gets the estimate 

(s-t) n a(s,t) n 



|U(s, t) :n E m \\ < CmCm+2 ■ ■ ■C m +2n-2- 



with the result that the series in (8) indeed converges to a unitary operator, for s — t small 
enough. Equation (7) then follows from (8) and allow us to extend the validity of the last 
conclusion — and, in turn, its own domain of validity. More detail on this is found in the 
important paper [18]. 

We can give an exact formula for the phase now. Use the standard identities 



|(«g,(t)) = 3|(logz(«)), | 



(logdet A(t)) = Tr A(t) 



-i dA(t) \ 
dt )' 



which give, from (5), 

d_ 

8X 



tt\T,t) = -\^(r u{T , t) 



d_ 



Therefore, 



Tr 



X = T 



t u(\,t) dr 



Tr 



U(\,t) 



d 



X = T 



Tu(X,T),Tu(r,t) 



dr. 



(The last expression is more symmetrical; the trace of this commutator is not zero, because 
it is taken in Fock space, whereupon, in view of the form of the scalar product (1), 



Tr 



A++ A+. 



A_ + A__ 
The phase of the scattering matrix is then 



TtA ++ + TtA__ .) 



6 = --9 
2 



Tr Tj 



■U(t, -oo) 



d_ 
dX 



A = T 



dr. 



(9) 



The analogous formula, with the same notation, for the boson case was first given, to the 
best of our knowledge, by Varilly and the author [19]; it differs only by a sign. Then, 
equivalent formulae both for the boson and fermion cases were found by Langmann [18]. 
The latter apply to charged fields, which are the ones considered in this paper. However, 
under the form (9) and with a suitable interpretation, the phase formula is applicable to 
Majorana fields, which are more general than charged fields [20] . 
Also note, before continuing, that 



d_ 

ds 



0(s,t) 



s=t 



d_ 

dX 



£(A,s,s) = 0. 



(10) 



A=s 



In other words, there is no contribution from the coincidence points of Tjj^x, t ) an d Tiy( T ,t)- 
This will prove to be the crucial remark. 

One has simply <9/<9A|, T[/(^ jT ) = iV dd(T~) in our present framework. Therefore, on 
calling T(r,t) := T u{T)t) , finally: 

6(8, t)=\J Tr (^+- (t)T_+(t, t) - T + _ (r, t)V-+(r)) dr, 

a rather elegant expression. Note that it differs from zero only at second order in pertur- 
bation theory. At that order, 

0*(s,t) = ±J Tr(V + .(T)U. +:1 (T,t)-U + . :1 (r,t)V. + (r))dT 
with an obvious notation. Since U : i(r, t) = —i V(r) dr, we set out to compute 
e :2 ( S ,t) = - l -J*Tr[V + _(r)^jy(X)dX^ - (jfV(A) dX^j VL+(r)] dr. 

Of course, (10) still applies at this approximation. The total phase 9-2 '■= O-^ipo, — oo) at 
this approximation is then 



0. 



OO POO 



oo J — oo 



8(t 1 -t 2 )(V + -(t 1 )V- + (t 2 )-V + -(t 2 )V- + (t 1 ))dt 1 dt2. 
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(11a) 



It should be clear that, at the same order of approximation, this is precisely 



~3< in 



/oo /-oo \ 
/ T[V(*i)V(* 2 )]d*id* 2 0iA (116) 
-oo J —oo I 



where T denotes the time-ordered product. 

4. Perturbative calculation of the phase in quantum electrodynamics 

In QED, from (2), or working like in the derivation of (3b), the integrals (11) are 
recast as 

e 2 f 

— 4 (xi)S + (xi — X2)4- { x 2)S~{x 2 — Xi)] d 4 xi d 4: X2. (12a) 

Here denote the Wightman "functions" . 
The first thing to remark is that, since 

S-(x)YS + (-x) - S+(x)YS-(-x) = S JP (x)YS + (-x) - S + (x)YSjp(-x), 

and Sjp has support inside the lightcone, then the integrand has support inside the light- 
cone. That allows one to substitute for 9{t\ — t 2 ) the covariant expression 9{{v{x\—X2))) =: 
x{xi — X2), where v is an arbitrary timelike vector, which can thus be varied at will, and 
a parenthesis has been used to denote the Minkowski product g^y^x 11 =: (yx) of two 
four- vectors y,x. Let then 

F» v (x) := tr[^S-{x)YS + {-x) - ^S + {x)Y S~ (-x)], 
F^{x) := X (x)F^(x). 

The second thing to remark is that the last expression — and hence (12a) — is only 
formal: is actually undefined as a product of distributions in view of the singularities 

of F^ v on the lightcone. Because the apparent trouble occurs at thejxiincidence points and 
since x(x)F Mi/ (x) makes sense for x 7^ 0, one can try to define xF^ as a distributional 
extension — or "regularization" in the terminology of [21] — of the latter. The scaling 
degree [5] or singular order [11] of the integral of the product of the two Wightman functions 
is 2; therefore, distinct extensions of this quantity will differ by linear combinations of the 
delta function at the origin and its derivatives up to order two — i.e., by polynomials in k 
of degree at most two in momentum space. The procedure is undoubtedly sound in the 
present case, as we know a priori the phase to be finite. Moreover, our framework will 
allow to select the "good" extension. 

It is indeed convenient to work in momentum space. There, 

0..2 = e - (2tt) 2 3 J F^(k)A^k)Mk) d 4 k, (126) 
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(27T) 



where we take into account that A{— k) = A(k), with the bar meaning here complex 
conjugation, because A(x) is real; and, formally, F^ v {k) = ( 2ir)- 1 / 2 X *F^(k) with * 
denoting ordinary convolution. We recall that, for timelike k and when choosing a frame 
in which k = (k°, 0), 

X(k) = = ( P 1 + i«6(*P)) . (13) 

To compute F^ u (k), one looks at the Fourier transform of tr[7 M S + (x)^ u S~ (x)]. By 
using the well known expressions of the Wightman functions in momentum space, this is 
expressed as 

— ^— r / tr[ 7 ^ + m)Y{4 - m)]6(p )6(q°)5(p 2 - m 2 ) 

7T) J 

x <% 2 - m 2 )5 4 (A; -p-q) d 4 q d 4 p =: -— l—T Vfl (k). (14) 

(2tt) 4 

Moreover, tr[(^ + m)7 M ((^ — m)7 l/ ] = A(p^q u + q^p v — ((pq) — m 2 )g flu ). Then one of the 
integrations in (14) is immediately disposed of, with the help of the 5 4 -function. The other 
is easily performed, with the help of the remaining 5-functions, again by choosing a frame 
in which k = (k , 0) so that F^ v can be regarded as a function of only one variable; we 
obtain 

T^(k) = f - <T J [k 2 (l + 7 (fc 2 ))(l - 2 1 (k 2 ))^ 2 9(l - 2 1 (k 2 ))0(k )} , 

where 7(/c 2 ) := 2m 2 /k 2 . All this is found in many books [11, 22, 23]. 

Thus we have been led formally to compute x * \ where F^ u {k) equals 

1 -(7^ [A; 2 (l+7(A; 2 ))(l-27(A; 2 )) 1 / 2 6»(l-27(A; 2 ))sign(A; )]. 



3(2tt) 3 V k 2 

That indeed behaves as a polynomial of degree two at high momentum transfer, confirming 
that the singular degree of F^ v is two. 

The correct (unique) recipe to regularize the imaginary part of x * is selected 
by prescribing that the result F^ v vanishes, together with derivatives up to order two, 
at zero momentum. This kills the delta function at the origin and its derivatives in con- 
figuration space, which otherwise would give a nonzero contribution to d/ds\ 6(s,t), 
contradicting (10). It is clear now, from the 5-function in (13), that 

UF^(k) = \F^{k), 

whereas the relation between the real and imaginary parts of F^ v is then given by a 
subtracted (at the origin) dispersion relation 



(fc ) 3 /" 00 UF^{Q 



( 3 ((-k°) 



d(. (15) 



The prescription that F^ v possess a zero of the third (indeed, fourth) order at k = 
can be independently justified by an argument that, although heuristic, we deem very 
strong in the present context. On invoking the Maxwell equations (in the Lorentz gauge) 
A^ik) = j^(k)/k 2 to conjure up the source j of the classical field, and on introducing 



G{k) := 1(1 + 7 (fc 2 ))(l - 2 1 (k 2 )) 1 ' 2 9(l - 2 7 (fc 2 )) sign(fc°), 



one gets 



y (2tt) 2 K J F^(k)A,(k)A v (k) d 4 k = -^-J (j(k)j(k))G(k) d 4 k. (16) 

The continuity equation (j(k)k) = has been employed to simplify the result. This simple 
expression exhibits only gauge-invariant variables. 

Now, we remark that (classical) gauge transformations are in general not imple- 
mentable in 1 + 3 dimensions. This is a very good indicator of the existence of ultraviolet 
divergences in the nonlinear theory, and indeed it was used by J. C. Varilly and the author 
in [3] to point out that QFT theories on noncommutative manifolds had to be ultraviolet 
divergent. On the other hand, in the linear theory selfinteraction is absent, so we would 
not expect ultraviolet divergences on physical grounds. That nonimplementability is the 
only source of spurious divergence difficulties. We therefore expect to be able to express 
the phase in terms of gauge-invariant variables, in a similar way to (16): 



0-.2 = ~ J (j(k)j(k))H(k) d% (17a) 

with H(k) is regular at k = 0; this is equivalent to F^ v having the aforementioned be- 
haviour at the origin in momentum space. 

Taking into account that G is odd, equation (15) leads immediately to a simple form 
for H: 

H(k)= i P r ( i7») 

TV J 4m 2 X(X-k z ) 

The restriction to timelike k can be removed by analytic continuation. Making the change 
of variable A =: 4m 2 / (1 — v 2 ), we get 



[ } ~ nk 2 J v 2 -l + 4m 2 /k 2 ' 



which is essentially the expression one finds in textbooks [24, pp. 249-252] for the vacuum 
polarization functional, after renormalization. The last integral can be easily carried out 
analytically, and it is then an instructive exercise to check that the function H is perfectly 
smooth at k = (at k 2 = 4m 2 , the onset of the absorptive part, H has a cusp). We shall 
not go into the details. 

Before rushing to the conclusions, a comment is in order: we have more or less treated 
the A(x) as test functions, guaranteeing implementability of the interpolating operators, 
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for the sake of the argument. However, it is clear that the final formula for the phase (9) 
is acceptable with only the milder requirement of the implementability of the scattering 
operator; consult [25] for a very efficient removal of technical conditions on the potentials, 
for this last purpose. 

5. Conclusions 

In this paper we have performed what amounts to an ab initio finite calculation of 
the "bubble" diagrams in linear quantum field theory. 

Now, at first significant order in QED this is essentially the same as the one-loop 
vacuum polarization or "photon self-energy" diagram (see in this respect [26, pp. 195-196]). 
The computation done here does not appear to have been pushed to that finish line before 
now, although the tools have been there since the seventies at least [27]. (Besides Varilly 
and the author [19,13], Langmann and Mickelsson [18,25] came close in the nineties.) The 
main point is that the "local causality condition" (10) selects the correct prescription 
among all the (finite) regularizations, with recourse to neither heuristic arguments [11] nor 
the extremely long and complicated "nonperturbative proof" in [28]. 

On the other hand, it will not have escaped the reader's attention that, in order to 
avoid pitfalls, we reorganize the calculation in the same way as [11]. The whole procedure 
is thus in the spirit of the Epstein-Glaser renormalization procedure, where there are Feyn- 
man graphs, but the Feynman rules do not necessarily apply — we avoided rewriting (12a) 
in terms of Feynman propagators. It is remarkable that, in the boson case, the quantum 
scattering matrix was found long ago by the Epstein-Glaser method by Bellissard [29] 
— without the phase. 

From formulae (16) and (17) one can easily verify a posteriori Bogoliubov's causality 
condition 



which was the starting point of Epstein and Glaser. 
Acknowledgments 

I am happy to thank Ricardo Estrada, Hector Figueroa, Steve Fulling, Carmelo P. 
Martin and Joe Varilly for discussions, and my students Esteban Araya and Pablo Blanco 
for help with drawing H. The warm hospitality of the departments of theoretical physics 
at the Complutense in Madrid and at Zaragoza, and the financial support of Sabaticos 
UCM 1999 are most gratefully acknowledged. 

References 

[1] A. Connes and D. Kreimer, "Renormalization in quantum field theory and the Rie- 
mann-Hilbert problem I: the Hopf algebra structure of graphs and the main theorem" , 



hep-th/9912092l IHES, 1999. 
[2] T. Filk, "Divergences in a field theory on quantum space", Phys. Lett. B 376 (1996), 





.o 

2' 



53-58. 



11 



J. C. Varilly and J. M. Gracia-Bondia, "On the ultraviolet behaviour of quantum 
fields over noncommutative manifolds", Int. J. Mod. Phys. A 14 (1999), 1305-1323. 

S. Minwalla, M. Van Raamsdonk and N. Seiberg, "Noncommutative perturbative dy- 
namics" , |hep-th/9912072| , Princeton, 1999. 

R. Brunetti and K. Fredenhagen, "Microlocal analysis and interacting quantum field 
theories: renormalization on physical backgrounds" , [math-ph/99U3U2^ , DESY, Ham- 
burg, 1999. 

D. R. Grigore, "On the uniqueness of the non-abelian gauge theories in the Epstein- 
Glaser approach to renormalization theory", [hep-th / 9806244| , Bucharest, 1998. 

T. Hurth and K. Skenderis, "The quantum Noether condition in terms of interacting 
fields" , |hep-th/981123Tj 1998. 

G. Pinter, "The action principle in Epstein-Glaser renormalization and renormaliza- 
tion of the S-matrix of $ 4 -theory" , |hep-th/99lHM3| , DESY, Hamburg, 1999. 

H. Epstein and V. Glaser, "The role of locality in perturbation theory", Ann. Inst. 
Henri Poincare A XIX (1973), 211-295. 

P. Blanchard and R. Seneor, "Green's functions for theories with massless particles 
(in perturbation theory)", Ann. Inst. Henri Poincare A XXIII (1975), 147-209. 

G. Scharf, Finite Quantum Electrodynamics: The Causal Approach, 2nd edition, Sprin- 
ger, Berlin, 1995. 

J. M. Gracia-Bondia, J. C. Varilly and H. Figueroa, Elements of Noncommutative 
Geometry, Birkhauser, Boston, 2000. 

J. M. Gracia-Bondia and J. C. Varilly, "QED in external fields from the spin repre- 
sentation", J. Math. Phys. 35 (1994), 3340-3367. 

D. Shale and W. F. Stinespring, "Spinor representations of infinite orthogonal groups" , 
J. Math. Mech. 14 (1965), 315-322. 

N. N. Bogoliubov and D. V. Shirkov, Introduction to the Theory of Quantized Fields, 
Wiley, New York, 1959. 

H. Grosse and E. Langmann, "A super-version of quasi-free second quantization: I. 
Charged particles", J. Math. Phys. 33 (1992), 1032-1046. 

J. T. Ottesen, Infinite Dimensional Groups and Algebras in Quantum Physics, Lecture 
Notes in Physics: Monographs m27, Springer, Berlin, 1995. 

E. Langmann, "Cocycles for boson and fermion Bogoliubov transformations" , J. Math. 
Phys. 35 (1994), 96-112. 

J. C. Varilly and J. M. Gracia-Bondia, "^-matrix from the metaplectic representa- 
tion", Mod. Phys. Lett. A 7 (1992), 659-667. 

F. Scheck, Electroweak and Strong Interactions , Springer, Heidelberg, 1996. 

R. Estrada, "Regularization of distributions", Int. J. Math. & Math. Sci. 21 (1998), 
625-636. 

[22] V. B. Berestetskii, E. M. Lifshitz and I. P. Pitaevskii, Quantum Electrodynamics (vol- 
ume 4 of the Landau and Lifshitz Course of Theoretical Physics), Pergamon Press, 
Oxford, 1980. 



12 



[23] C. Itzykson and J.-B. Zuber, Quantum Field Theory, McGraw-Hill, New York, 1980. 

[24] W. Greiner and J. Reinhardt, Quantum Electrodynamics, Springer, Berlin, 1992. 

[25] E. Langmann and J. Mickelsson, "Scattering matrix in external field problems", J. 
Math. Phys. 37 (1996), 3933-3953. 

[26] J. M. Jauch and F. Rohrlich, The Theory of Photons and Electrons, Springer, Berlin, 
1976. 

[27] S. N. M. Ruijsenaars, "Charged particles in external fields. II. The quantized Dirac 
and Klein-Gordon theories", Commun. Math. Phys. 52 (1977), 267-294. 

[28] G. Scharf and W. F. Wreszinski, "The causal phase in Quantum Electrodynamics", 
Nuovo Cim. A 93 (1986), 1-27. 

[29] J. Bellissard, "Quantized fields in interaction with external fields I. Exact solutions 
and perturbative expansions", Commun. Math. Phys. 41 (1975), 235-266. 



13 



